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Abstract. Following the Nambu mechanics framework we demonstrate that the 
non dissipative part of the Lorenz system can be generated by the intersection 
of two quadratic surfaces that form a doublet under the group »S'L(2,R). All 
manifolds are classified into four distinct classes; parabolic, elliptical, cylindrical 
and hyperbolic. The Lorenz attractor is localized by a specific infinite set of one 
parameter family of these surfaces. The different classes correspond to different 
physical systems. The Lorenz system is identified as a charged rigid body in a 
uniform magnetic field with external torque and this system is generalized to give 
new strange attractors. 
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1. Introduction 

We discuss the geometric structure of the integrable part of a chaotic three dimensional 
dynamical system. Our approach is generic, however we concentrate on the Lorenz 
model p!^. The analysis is performed in the Nambu mechanics [2 formalism. We find 
that important information regarding the boundary of the strange attractor can be 
extracted from this geometric structure. In addition, the symmetry that underlies the 
phase space geometry leads to mappings between different dynamical systems. Let us 
review some basic concepts. 

In [2 Y. Nambu introduced a formalism of classical mechanics to three 
dimensional phase spaces. In his words "the Liouville theorem is taken as the guiding 
principle". Hence, the formalism applies to systems that preserve the phase space 
volume, i.e. \/v = 0. The motion of points x = (x, z) G in phase space is 
determined by two functions of {x^y^z), called Nambu Hamiltonians Hi and i^2, by 
the equations: 

f=Vi^iXVi^2 (1) 

By definition ([l]) the fiow is volume preserving, since Vv = 0^ v = x. In general any 
function z) is evolving according to 

F = VF • {\/Hi X \/H2) = e^jkd^FdjHidkH2 (2) 

where Sijk is the Levi-Civita tensor. The right-hand side defines the Nambu three 
bracket: 

{F, HuH2} = e,jkd^FdjHidkH2 (3) 

The Nambu Hamiltonians are conserved quantities as one can see from ([2| . The phase 
space orbit can therefore be realized as the intersection of the surfaces Hi = Hi{t = 0) 
and H2 = H2{t = 0). 

In [3, 4 the Nambu bracket is realized as the Poisson bracket, denoted as {•, -ji/s' 
induced on the manifold H2 = H2{t = 0), which is now upgraded to a phase space. 
We will call this bracket Nambu- Poisson bracket and is defined as in (|3|. It is 
antisymmetric and satisfies the Jacobi identity. The Nambu-Poisson bracket of one 
dynamical variable Xi and Hi is: 

{xi,Hi}H2 = eqjkdqXidjHidkH2 = eijkdjHidkH2 = VHi x VH2 
The dynamical variables form an algebra on H2 

{xi, Xj}H2 = £qpkdqXidpXjdkH2 = eijkdkH2 (4) 
The evolution of the system is given by 

Xi = {xi,Hi}H2 (5) 

The surface H2{x^y^z) = const. = H2{t = 0) can in some cases be regarded as the 
group manifold of a Lie group, for example if = + + this specific 

example H2 is just the Casimir of 5*0 (3) , so it is reasonable to call it C and we will 
do so often in this work, calling Hi as just H as well (in some cases we will stick 
to Hi^ H2 notation). Now, C defines an algebra of generalized coordinates x^y^z by 
equation Q: 

{x,y}c = z 
{y,z}c = X 
{z,x}c y 
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The evolution lies on surface C = C{t = 0) = const., with dynamics determined by 
H , or equivalently one can think of the operator {•, H}c as generating the flow by 
repeated operations of a group element of the group with Casimir C (which group 
element is determined by H). Therefore, in this case, the surface on which the system 
is constrained is a group manifold. 

There is a freedom [2_ in determining Hi, H2 since any transformation 
{H{,H2) = {H[{Hi,H-2),H2{Hi,H2)) with Jacobian determinant equal to one leaves 
the equations of motion invariant: 

VH[ X VH'2 = eijkdjH[dkH'2 



- '^'^ \dIhdH2 ~dH2dIh) ^^^'^'^^^ 

We will call this freedom a gauge freedom. In the case of linear transformations 
the gauge group is SL{2^R) and therefore the Nambu Hamiltonians form a doublet 
h = {Hi^H2) that transforms under 6'L(2,R). 

In this work we will use a decomposition of the Lorenz system ^ , which has three 
phase space dimensions, to a conservative (non-dissipative) part and a dissipative part 
[5]. In Hamiltonian mechanics (even phase space dimensions) it is common to add to 
some known Hamiltonian physical system, forced and dissipative terms. However, it 
is rarely the case that given the equations of motion of a forced, dissipative system 
{Vv 7^ 0), one is asked to guess a decomposition in a Hamiltonian (non-dissipative) 
part and a forced, dissipative part. This decomposition is of course not unique, but 
in Hamiltonian mechanics it is often simple and natural to guess 'the most physical' 
decomposition. For example, suppose one is given the system 

Xi = X2 

±2 = —uoxi — bX2 

which is a dissipative system (Vv = —b). It may be decomposed to a non-dissipative 
part vnd with \/vnd = and a dissipative part vd- 

V vnd + vd 

One may choose the parts vnd = (0, — cjxi) and vd = (x2,— 6x2) or even vnd = 
{x2 +X2,— cjxi) and vjj = (— X2,— 6x2). Mathematically they are nice equivalent 
choices but as far as physics is concerned 'the most physical' choice would have been 
^ND = (^2, —ooxi) to get the harmonic oscillator for the non-dissipative part and to get 
the friction force F = —hx {velocity) for the remaining dissipative part vd = (0, —6x2). 

The situation is analogous in three phase space dimensions using Nambu 
mechanics, even though one is not guided by intuition to find a physical decomposition. 
The first time a decomposition in Nambu mechanics was applied was in [5 (to the 
Lorenz system). In [6j has been set a general framework of the dissipative Nambu 
mechanics with more examples. The decomposition is chosen in [6 so as \/vnd = 
and \/ X Vd = 0. It is not unique, however one is restricted by convenience according 
to motivations. The velocity field can always, at least locally, (see [6 and references 
therein) be written as 

v = \/x A^\/D , with Vnd =S/ x A and vd 
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for some vector field A and some scalar field D. The field A has a gauge degree of 
freedom a{x) 

A^ A^Va 

Any vector field can be written A = Va + [8] for some scalar functions a, /3, 7. 
Working with specific Hamiltonians Hi and H2 means that you have chosen the gauge 
A = H1VH2 for P = Hi and j = H2. This gauge is known as the 'Clebsch-Monge' 
gauge and the scalar functions a, /3, 7 are usually called 'Monge's potentials' (see also 
[71 [9j for more details on the application of Nambu mechanics on non-Hamiltonian 
chaos and the quantum, non-commutative approach). 

We are interested, and will focus in this work, in the gauge freedom 5^(2, R) 
of Hi and H2. The application is on the Lorenz system, due to its unique 
significance regarding chaotic dynamics and the interesting physical properties of its 
non-dissipative part in the specific decomposition we will use. We classify the 'Nambu 
doublets', i.e. the pair of Hamiltonians, of the Lorenz non-dissipative part with 
criterion the isomorphism of the algebras. We find four different classes of doublets in 
section [2] that are identified with four different systems in section [4| which are therefore 
dynamically equivalent. Three of them are very familiar physical systems, namely the 
simple pendulum, the Duffing oscillator and a charged dielectric rigid body inside a 
homogeneous magnetic field, with the fourth being an 50(2, 1) system. We may say 
that these four systems are 'dual' to each other. 

We stress that Lorenz defines his system as a forced, dissipative system [1 . He 
starts from an abstract, conservative system with some conservative quantity Q which 
can be a compact 'shell like^ surface. Then he makes the hypothesis that '[•••], 
whenever Q equals or exceeds some fixed value Qi, the dissipation acts to diminish 
Q more rapidly then the forcing can increase Q, then {—dQ/dt) has a positive lower 
bound where Q > Qi, and each trajectory must ultimately become trapped in the region 
where Q < Qi.' quoted straight from the famous 1963 paper [1^. This reasoning 
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(a) The phase space trajectory of a damped 
harmonic oscihator. With red the non- 
dissipative trajectory (harmonic oscillator) 
corresponding to the initial conditions. 
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(b) The Lorenz system for large initial 
(Xi, X2, Xs) projected on Xi — X2 plane. 
The spiral is evolving towards decreasing 
Xs. With red the non-dissipative trajectory 
corresponding to the initial conditions. 



Figure 1. The damped harmonic oscillator and the Lorenz system. In the Lorenz 
system, despite of being a dissipative system (Vv < 0), there is a constant driven 
term that appears in specific coordinates. For large Xi ^s, however, the Lorenz 
system behaves as a purely dissipative (not driven) system. 
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implies that the conserved quantities of the corresponding conservative system, which 
in our decomposition are the Nambu function^ is most natural to be used for finding 
trapping surfaces and for the localization of the attractor in general. This is what we 
do in section [3] using the Nambu functions of section [2] 

In section |4.3| we find that the Lorenz system is equivalent to a physical system 
that presents this type of behaviour, quoted earlier, for the forced, damped terms. This 
system is the SO (3) system we mentioned (charged rigid body in a magnetic field), 
but now with external torques. The generalized coordinates are the three components 
of the angular momentum L. The conservative function "Q" of Lorenz is just the 
measure of angular momentum Q = \L\ and its conservation for the conservative part 
expresses the conservation of angular momentum. The damped terms are torques 
proportional to angular momentum, like air resistance. There is a constant torque 
along Z-direction that acts as a driven (forced) term when < 0. For large L the 
friction terms dominate over the constant torque (see Figure [T]). The attractor in this 
coordinates lies in Lz < 0. Since there exists a constant driven term for L^, the 
angular momentum could not ever become constantly zero. 



This system ( 40 ) is generalized for arbitrary direction of the magnetic field to give 



a new dynamical system. This general system (44) presents new strange attractors 
given in Figures [TT] [12] 



2. Lorenz non dissipative part and 6'L(2,R) classification 

The Lorenz system [1 , defined by: 

X = aY -aX 

Y = -XZ + rX -Y (7) 

Z = XY -bZ 

where r, a and b are positive constant^ is a dissipative system, since \/v = 
— (cr + 1 + 6) < where v = (X, F, Z). In [5 the following decomposition is proposed: 

VND = {(tY, -XZ + rX, XY) , vd = (-aX, -F, -bZ) (8) 

with the full system given hy v = vnd -\- ^D- In the same work two Nambu 
Hamiltonians for the non-dissipative part vnd are found 

H =^Y^ ^^Z^ -rZ , C = -^X^ + aZ (9) 

by integration of the equations: 

dY _r-Z dZ _X 
lz ~ Y ' Zx ~ 7 



(10) 



To decouple the parameter dependence between the two parts (non-dissipative and 
dissipative) and for reasons of convenience for our purposes we introduce the scaling 
transformation: 



X = Xi , r = V^/aX2 , Z = ^/r/aXs , p = (11) 

J We sometimes call the Nambu Hamiltonians of the non-dissipative system, just 'Nambu functions'. 
§ In every figure of the Lorenz system the free parameters are assumed to be r = 28, cr = 10 and 
6 = 8/3, unless stated differently. 
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and the Lorenz system ^ becomes: 

Xi - pX2 -aXi 

X2 = -X1X3 + pXi -X2 (12) 
Xs = X1X2 —bXs 
We will work independently on the non dissipative part: 
Xi pX2 

X2 = -X1X3 + pXi (13) 
Xs = X1X2 

The Nambu Hamiltonians ([9| become in the new rescaled variables: 

H = \xl + ]^xl-pX:i,C=-\xl+ pXs (14) 
We will call this pair of functions a 'Nambu doublet' 

h=(^] (15) 



since it transforms as a whole under a restricted class of transformations. These are 
all transformations with Jacobian determinant equal to one, as we saw in section [T] 
We will restrict ourselves in linear transformations so as to work only with quadratic 
functions. In this case, the transformation group is the group of real 2x2 matrices 
with determinant one, namely the 5^(2, R). For example let generate a different 
doublet by action of an S'L(2,R) matrix: 

We see that the role of H and C can be interchanged, as long as you change the sign 
of one of the functions, since we could have used 

-1 



^~ ' 1 
as wellTll Another transformation could be: 



Co I 1 I [ C 



that gives: 



^0 J V "2^1 + 2^2 + 2^3 



(17) 



Recall from the Introduction, that the Nambu Hamiltonians i^, C of a volume 
preserving system define surfaces in phase space because they are conserved. Once 
initial conditions are given, the solution of the system should lie on both surfaces 
defined by 

H{XuX2,Xs) \_f i/(Xi(0),X2(0),X3(0)) 



c(Xi,X2,X3) ; V c(Xi(o),X2(o),X3(o)) ; ^^^^ 

Therefore the intersection of the two Hamiltonians defines the phase space trajectory 
(see Figure [2|. Thus, given a Nambu doublet and initial conditions, it is not only 
given the equations of motion (velocity field) but the solution of the system, as well! 

II It is apparent that the system is invariant in these cases since Xi = VH x VC = (— VC) x VH = 
VC X i-VH) 
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(a) 



(b) 



(c) 






(d) 



(e) 



(f) 



Figure 2. An orbit, in red, of the conservative system ( |13| > as defined by the 
intersection of various Nambu surfaces (same orbit at each case). 



In equation (17) is defined the Nambu Hamiltonian function Co that defines a 



hyperboloid, once initial conditions are given. This is a topologically different surface 
from the cylinder and parabolic surface of equation (14). An /SL(2,R) transformation 



led to a new type of surface. This is reflected to the algebra of the Nambu bracket, 
since: 

{Xi,X2}c = P 

{X2,X^}c = -X^ (19) 
{X3,Xi}c = 



while: 



^3 
X2 



(20) 



{Xi,X2}co 
{X2,X3}co 
{X3,Xi}co 

The later algebra is the 50(2, 1) algebraj^ We could get an isomorphic to 5*0(2, 1) 
algebra by applying to (17) transformations of the form 

(21) 



A: 



1/C 



for some C G K. These algebras would correspond to other hyperboloids: 



C 



1 1 

'2C 



y U sing the algebra ( |20t, the properties of the Poisson bracket and Hq one can show that the doublet 
\17\ gives the system^Ts} : 



Xi = {Xi,Ho}co = {^1, |^?}co +{^i,-P^3}co 

X2 = {X2,Ho}co = {^2, |^?}Co + {^2, -pX3}co 
X3 = {X3,Ho}cq = {X3, + {-^3, -P^sjco = ^1^2 



-9X2 
■ —X1X3 - 



-pXi 
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One may ask what other types of quadratic surfaces one may get. We answer this 
question and classify the doublets according to isomorphisms of the algebras. 
Let an arbitrary 5^(2, R) matrix 



A 



5 



aS — /3j ■■ 



act on the doublet (17). Then we get a general expression for the possible quadratic 
doublets 

' 5(Xi,X2,X3;a,/3) " 
S{X,,X2,Xs'n,S) 

with 



h- 



ya,P,-f,5 eR: a5 - = 1 (22) 



5(Xi,X2,X3;a,/3) = (a- 



apXs (23) 



Of course it does not matter with what doublet one may begin, since one can 



always with a re-parametrization end up with equation (23). The expression for 
6'(Xi, X2, X3; 7, (5) is identical with 23 and describes exactly the same surfaces. The 
S{a^P) and S{'y^S) are different only when considered as part of a doublet, because 
it is then imposed the restriction aS — = 1. 



(A) 



A < 1 

A > 1 

A = 1 

A =b(X) 



(A-l)|X? + ix| + |(X3-Ap)2 

a set of hyperboloids centred at X| = Xp 

a set of ellipsoids centred at X| = Xp 

a cylinder centred at X| = p 

a parabolic cylinder 



Table 1. The Nambu surfaces 



The possible types of quadratic surfaces are given by equation (23). These 



determine the possible Nambu algebras and the geometry of the phase-space through 
the Nambu bracket: 



Let A = 



{Xi,Xj}s 
and 



(A) 



(A 



(24) 



where is an abbreviation for C(Xi,X2,X3; A). Then, up to a constant shift it is: 

S{X^,X2,X:,■,a,0) = pC(X) 

Rescaling of a Nambu function: H[ = PHi leaves the corresponding Nambu surface 
invariant since: H[ = H[{0) <^ PHi = PHi{0) 4^ Hi = i^i(O) so that the equations 
H[ = H[{0) and Hi = Hi{0) define exactly the same surface in the phase space. 
Therefore equation (24) describes all possible quadratic surfaces, which are listed in 



+ Note that there does not exist another different distinct set of quadratic doublets. If it existed 
would mean that members of one set are not related with members of the other set by transformations 
with Jacobian determinant one, which is not allowed as can be realized by (|6|. In addition quadratic 
doublets can only be related with linear SL(2,M.) transformations with each other, because any 
non-linear transformation would raise the rank of the functions. 
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(c) Cylindrical (d) Elliptical 

Figure 3. The four different classes of the Nambu surfaces (invariant manifolds) 
of the non dissipative part ( |13[ ) of the Lorenz system. An orbit, in red, of the 
system lies on any of these surfaces. 



Table [T] and can be seen in Figure [3] From the algebraic point of view the rescaling 
leads to isomorphic algebras. We see that there are infinite surfaces, of four types, 
when classified with respect to the isomorphisms of the corresponding algebras. The 
surface given by 0) (that is C of equation (14)) corresponds to A ^ ±oo. 

Let list the possible algebras. The hyperboloid centred at the origin Cq gives the 
50(2, 1) algebra and any other hyperboloid gives an algebra isomorphic to 5*0(2, 1). 

The sphere O2 = ^Xf + + ^(Xs - 2^)^ centred at = 2p, after the 
transformation X3 ^ X3 — 2/), gives the 50(3) algebra {X^,Xj}c2 = ^ijk^k , and 
any ellipsoid, an algebra isomorphic to 50(3). 

The cylinder d = |X| + |(X3 - p)^ gives the SE{2) algebra after the 
transformation X3 X3 — p. 

The algebra obtained by the parabolic cylinder 0(oo) call it O = + pX^, 

can be realized as the Heisenberg algebra {Xi, X2}c = p of two independent variables 
Xi,X2 on the curved plane O: {X2,X3}c' = — Xi , {X3,Xi}c' = 0. 

With these four types we classify the possible Nambu doublets as in Table |2] and 
an orbit defined by six different doublets can be seen in Figure [2I In Table [2] one of 
the four possible types of algebras is chosen with a specific C at each case and the 
corresponding H is found by the condition a5 — {3'^ = 1. 
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H = {a- l)^Xf + + |X| - apX3 , Va G M 
{Xi,X2}c = P 

{X2, Xsjc" = —Xi 
{X3,Xi}c = 

if = (l-/3)|x2+/3|x|+/3|x|-pX3 , V/3GM 
{Xi,X2}cq = X^ 
50(2,1) {X2,X3}co =-Xi 

H = \XI + + - + 1)X3 , G R 

{Xi,X2}ci = ^3 - P 
5£;(2) {X2,X3}ci=0 

{X3, Xijc^ = X2 
if = (/3 + + ^\xl + - p(2/5 + 1)X3 , V/3 G : 

^0(3) {X2,X3}C2 =^1 
{X3, Xi}c2 — ^'2 

Table 2. The Nambu doublets Classification 



3. Localization of the Lorenz Attractor 



When the dissipative part vd = (— crXi, — X2, —bXs) is added to the system (13), the 
Lorenz attractor is formed. This is a global attractor, that is this subset of phase space 
that every trajectory of any initial conditions approaches for t ^ 00. By definition the 
attractor should lie in a bounded region and determining this region is what we call 
localization. The possible geometries of the boundary of the Lorenz attractor have 
been discussed by Lorenz [1] [10] , Sparrow [11] and others [121 US] • We will not find 
new surfaces, not mentioned in bibliography already. Our scope is only to prove that 
the Nambu functions found in previous section, i.e. the quadratic invariant manifolds 
of the non-dissipative part, define proper localization surfaces. 



Let us change the notation and use just 5* for 2C(a) of equation (24). 

S{XuX2,Xs) = (A - l)X^^Xl + (X3 - Xpf (25) 

When Xi, X2, Xs are subject to the Lorenz evolution (with the dissipative part 
added), 5* is not conserved but is changing with time: 

S = ^b{Xpf - 2 (^a(A - + Xl + b{Xs - \xpf^ (26) 

Assume for the moment A > 1. Suppose there exists some real constant Kmax for 
which S can be positive only inside the ellipsoid 5* = Kmax and nowhere else (see 
Appendix C of Sparrow [11 and Doering and Gibbon [l2l, as well). Then for any point 
outside this ellipsoid, 5* < and therefore 5* will be decreasing until it gets S < Kmax 
i.e. the trajectory gets inside the ellipsoid. Of course, if the trajectory starts inside 
this ellipsoid, it will remain there, since on the outside is everywhere 5* < 0. If such 
Kmax exists, it is the solution of the constrained maximization problem: 

Kmax=max{S{X^,X2,X^)'.S>0} ,A>1 (27) 



This is a typical problem to solve with the Lagrange multipliers method (Appendix 



ffl. This condition defines a compact attracting surface (in fact every other ellipsoid 



Phase Space Geometry and Chaotic Attractors 



11 




Figure 4. The repelling two-sheeted hyperboloid A = 0.8 and the Lorenz 
attractor. 



S = K with K > Kjnax is an attracting surface, since 5* < everywhere on the outside 
of ah of these). The solution of (27) is given in (30) where K^ax = R^- 



We argue here that similar reasoning can be applied for non-positive-definite 
functions (non-compact surfaces). For A < 1 the surfaces are hyperboloids which are 
non compact. In this case the equation 6* = for some constant K defines a one 
sheeted hyperboloid if K > and a two sheeted hyperboloid if K < 0. For simplicity, 
with no loss of generality, assume that K < 0. Then the condition S > K defines 
the region outside the 'lobes' of the two sheeted hyperboloid (the attractor lies in this 
region in Figure [4|, while S < K defines the region inside the lobes. As previously, 
assume there exists some constant value Kmin such that S can be negative only outside 
the lobes (region S > Kmin) and nowhere else. Then any trajectory starting from 
inside the lobes will be expelled to the outside, since 5* > in the inside and therefore 
S is increasing until it gets S > K^in i-e. the trajectory gets outside the lobes. Of 
course, if it starts outside the lobes, it is impossible to get inside. We call such a kind 
of surface, a repelling surface. The same reasoning can be applied to any non-compact 
surface. The problem is now a minimization problenj^ 

Krmn=min{S{XuX2,Xs):S<0} ,A<1 (28) 

The solution is given in (30) where Kmin = —Rr- In the Appendix A is given 
analytically the Lagrange's multipliers method. 

Another way to understand localization is to think of a fixed surface and imagine 
the fiow at each point crossing it (Giacomini and Neukirch [13j). For a localization 
surface the fiow should cross every point of the surface towards the same direction 
(semi-permeable surface) inwards for the attracting (compact) surfaces and 

outwards for the repelling (non-compact) ones. Or else the trajectory could pass 
through the surface in one direction at some point and get through the opposite 
direction at some other point of the surface. Thus, the scalar product of the vector 
normal to the surface \/S at each point and the vector tangent to the fiow v should 
have a constant sign at each point. That is: 

S = (V^) • X (29) 

* The signs is a matter of convention, since we could have used a function S = —S and turn the 
problem to a maximization problem. 
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Figure 5. The homoclinic orbits of the conservative system \13\ with red as 
defined by the repehing cone and parabohc cyhnder. When the dissipative part 
is added, the flow is forced towards the Lorenz attractor. 



has constant sign at every point on a localization surface. This is proven in the 
Appendix B for all repelling surfaces of Equation (31). 



Let summarize our results. We find the following attracting and repelling surfaces: 

S{Xi,X2,Xs) = K<=>{X- l)Xf + X| + (Xs - Ap)2 = K 
>Rl for A > 1 



\JK : K { 



where: Ra 




for {A = 0, -oo} and {A = (-oo, 0), X3 > 0} 

for A < 1 

for 6 < 2 



(30) 



for 6 > 2 



, Rj. — \p , Rj. — \p 



The attracting cases hold for a > 1, the repelling cases hold for a > 1 and cr > h except 
for the parabolic surface A ^ — 00 that is —X^ + 2pX^ = which holds for <j > b/2. 
The surfaces that are closer to the attractor and provide the optimum localization are 



the extreme cases of (30) that correspond to the equalities: 

-Rr 



S — Rfj^ , S — R^ , 



S = 



(31) 



This means that the strange attractor lies in the region A that is defined as: 

A = {(Xi, X2, X3) :S<RldindS>Rl and S > -R^} (32) 

The following surfaces of (30) can be found in Doering and Gibbon flT: the 
repelling parabolic cylinder — Xf + 2pX3 = (A —00), the repelling cone 
-Xf + X| + X| = (A = 0), the attracting cylinder X| + (X3 - p)^ = (A = 1) 
and the attracting sphere X^ + X| + (X3 - 2/?)^ = (2p)^ (A = 2). The two-sheeted 
hyperboloids (A < 1) (see Figure |4| were for the first time found by Giacomini and 
Neukirch [13j. In Figure [s] we see that the repelling cone and the repelling parabolic 
surface define the homoclinic orbits of the non-dissipative system. For b < 2 these 
homoclinic orbits define also the attracting surfaces (Figure [6|. 
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\ -40 



Figure 6. For b < 2 the strange attractor is localized by the corresponding 
surfaces of the homoclinic orbit. Here for 6 = 2 are shown just three of these 
surfaces. The strange attractor is fairly seen in blue in the interior and the 
homoclinic orbit is shown in red. 



4. Correspondence with physical systems 

The different Nambu doublets can be interpreted as different physical systems. The 
one Nambu Hamiltonian is regarded to be the energy of the system and the other 
has different interpretations depending on the case (a constraint of the motion or an 
integral of motion) . This framework is providing a way to identify completely different 
physical systems that have the same dynamics. From the mathematical point of view 
it is just one system in different representations. In our case these are the Heisenberg, 
SE{2)^ S0{3) and 50(2, 1) formulations that are linked with each other. 

4.1. The unforced, undamped Duffing oscillator 

Let us work with the first doublet (i^, C) of Table |2j We reduce the system on the 
parabolic cylinder C = — + pX^ and the Hamiltonian becomes: 

H = \xl - ^(P' - C)Xl + + const. 

We see that VA the Hamiltonian is the same up to a constant shift. Let p be absorbed 
in the Hamiltonian, then 

^=\pn-\-y-C)Xl^l-Xt (33) 

Our dynamical variables are now the conjugate canonical variables (Xi,X2) and the 
Nambu bracket reduces to Poisson bracket on a 2D phase space (Figure [7|. 

{Xi,X2}c = l 

The system ([33| describes a particle moving along X\ with momentum under the 
influence of the potential: V[Xx) = — ^{p'^ — C)Xf -\- -^Xf where C = const, depends 
on the initial conditions. The equilibria depend on the value of C. They are X = (0, 0) 
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Figure 7. The phase space portrait of the system ( |13| red uced on the parabohc 
cyhnder C. Choosing initial conditions for the system \13\ is equivalent to 
choosing not only initial conditions but also a potential for the system ( |33| . The 
value of the potential is determined by the level of the surface C. 




Figure 8. The phase space of system \33\ . It depends on the initial conditions 
of system \13\ and specifically on Xs{t = 0) through the value of the function C. 



for C > and 



X 



{ (0,0) 



C),0) 



for C <p^ 



(V2(P^-C),0) 



This system is the unforced, undamped Duffing oscillator x = —ax—^x^ for a = C—p^ 
and P = \- For C > p^ the potential is a single well while for C < p^ the potential is a 
double well and the symmetry of the vacuum state is 'spontaneous broken' (Figure |8|. 

This simple analogue of a mechanism of 'spontaneous symmetry breaking' gives an 
intuitive picture of the formation of the attractor. Equation C — p^ defines a critical 
parabolic surface that separates phase space in the upper region of one vacuum state 
and the lower region of two vacuum states (Figure [9(a) ). Suppose the initial conditions 
are above that surface and dissipation is "turned on" . The system will flow towards 
the origin until it passes the critical surface when it enters the broken symmetry region 
and has to choose between different vacuum states. It may oscillate around a while but 
eventually it is forced outside the critical surface again. The procedure is constantly 
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(a) The 3D phase space of the non 
dissipative part ( |13| ) and the critical 
paraboHc surface. 



(b) The Lorenz system for initial 
conditions X = (0, 20, 40) and the 
critical parabolic surface. 



Figure 9. The resemblance with a 'spontaneous symmetry breaking mechanism' 
gives an intuitive explanation for the formation of the attractor. 



repeated and the attractor is formed (see Figure [9|. 

This is a quahtative picture. However, when the dissipative part is added, 
the potential becomes another independent variable of the system. With the 
transformation X2 — X2 — ^Xi the full Lorenz system becomes: 

Xi = pX2 

h = -j-Xl ^\{p^-cj- C)X, - (a + 1)X2 
C ={(j-\)Xi-hC 

We see that a locally time dependent potential cannot be defined. The oscillator 
develops a friction term and we obtain the Takeyama [I4 memory term in the 
potential: C{t) = C{0)e-^* + (a - |)e-^^ d^e^^Xf{^). 



4.2. Simple pendulum 

Now, let work on the Nambu doublet of Table |2] corresponding to the cylinder 
Ci = + ^(-^3 — p)'^' Let apply the transformation X3 = X3 — p. The non 

dissipative part (13) becomes: 

VND = (pX2, -X1X3, X1X2) (34) 

and the full Lorenz system becomes: 

V = {pX2 - aXi, -X1X3 - X2, X1X2 - bXs - bp) (35) 

The Lorenz system is just translated along X3 axis with the equilibrium of the origin 
translated at X3 = —p. Let study independently the volume preserving part (34). 
The cylinder becomes: Ci = ^Xl + ^Xf and we choose the Hamiltonian for /3 = 0: 

H = \xl - pX3 (36) 

Note the analogy of this Hamiltonian to the one of a simple pendulum 
\pQ/{ml'^) — mglcosO. We perform the following transformations: 

Xi = — X2 = RsinO^ X3 = RcosO 
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Figure 10. The phase space portrait of the sy stem ( |13| ) reduced on the cyhnder 
Ci. Choosing initial conditions for the system ( |13| is equivalent to choosing both 
initial conditions and the length for a simple pendulum. The length is determined 
by the radius of the cylinder. 



where our new variables are (O^po) since: R = -y/2C(i) = const. We see that Ci is 
just expressing the constraint that the simple pendulum has constant length. In fact 
we have parametrized Ci with (O^pe): Xi = XAO^^po). We have: 6 = tan~^^ and 
therefore 

dO _ dO Xs dO X2 
^1 



dXi ' 8X2 2Ci dXs 2Ci 

Using these derivatives is straightforward to calculate the following Nambu-Poisson 
brackets: 

OH OH 
{0,pe}ci 1, {0,H}c^ -g^, {pe,H}c^ ~~00 

Hence, the canonical equations are retrieved with: 

H = -pQ — pR cos 



Thus, Lorenz volume preserving motion (34) reduces to pendulum motion [|| on level 
surfaces of Ci (Figure 10). 

When the dissipative part is added the equations of motion (full Lorenz system) 
become: 

RO = Rpe-{l- b)Rs'm 6 cos 6- bp sin 6 

pe = —pRs'mO — ape (37) 
R = -bpcosO - R{sm^ + b cos^ 0) 

4-3. Charged rigid body in a uniform magnetic field 

Let us work on the Nambu doublet of C2 = \Xl + ^X^ + ^(-^3 — 2p)^ with the 
transformations 

Ll=Xi, L2=^2, ^3 = ^3-2/9 (38) 

jl An analogous correspondence between the simple pendulum and the free rigid body can be found 
in [16] 
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Figure 11. The phase-space of system for h = 1/8, h = 1/3, I3 = 1/2, 
pi = 0, p2 = 22, p3 = 30, cr = 10, 6 = 1 and initial conditions Lq = (3, —4, —10). 



The non dissipative part becomes: 

vnd (p^2, -LiLs 
and the fuh Lorenz system becomes: 
Li pL2 

L2 = — I/1L3 — pLi 
Ls = L1L2 



■ pLi,LiL2) 



-aLi 



(39) 



(40) 



2bp 



The Lorenz system is just translated along Xs axis with the equilibrium of the origin 
translated at L3 = —2p. The Nambu function C2 becomes: 



Co 



2 ' 2 ^ 2 

giving the algebra of SO {3). Let work with the corresponding Hamiltonian for (3 



H = Ll 



1 



(41) 
= 1: 

(42) 



We identify Li as the angular momenta and C2 is expressing the conservation of \L\. 
This system can be realized as a uniformly charged dielectric rigid body with charge 
Q in a uniform magnetic field B = (0, 0, B) with moments of inertia along the primary 
axis /i = 1/2 and I2 = I3 = I and p = QB. The dissipative part is just the external 
torque: a damped term proportional to angular momentum for each component and a 
constant term —2bp in L3 component that acts as a driven torque for L3 < 0. For an 
arbitrary rigid body with any moments of inertia and any direction of the magnetic 
field the Hamiltonian is: 



where 



H — Hfree + ^int 
Hfree = 277^1 ^ 

Hint = —QB • L 



1 r2 



1 r2 
2T^^3 



(43) 
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Figure 12. The phase-space of system ( [ii} for h = 1/20, I2 = 1/10, I3 = 1/5, 
pi = 10, P2 = P3 = 25, cr = 10, 6 = 1 and initial conditions Lq = ( — 1, 1, —1). 



Assuming p = QB and adding the external torque N^xt = (— cr^i, —^2, —^^3 — 26p) 
we get the equations of motion for the general system: 









= {L,,H}c,^ 


- Next 




^1 = 1 


^ 1 


1 ^ 
h 1 


-\-psL2 


- P2L3 


—aLi 




^ 1 


1 \ 
h J 


L1L3 +P1L3 


- P3L1 




^ = ( 


^ 1 


1 ^ 
/2 ; 


L1L2 -\-p2L1 


- P1L2 


-bLs - 2bp3 



(44) 



which can also be written as L-\-ujx L-\-QB x L = Next- System (44) presents strange 



attractors as can be seen in Figure and Figure 12 



4-4- SO (2^ 1) formulation 

The 5'0(2, 1) case can be regarded as a 'hyperbolic' analogue of the rigid body 50(3) 
case. This is mainly of mathematical interest. One gets this hyperbolic analogue with 
the Nambu doublet (/3 = 2 for example) 

Co = -\xl + \xl + \xhH = - \xl + Xl + Xl - pX3 (45) 

Depending on the initial conditions the hyperboloid Co may be two-sheeted which 
corresponds to the hyperbolic space or one-sheeted, which corresponds to anti-de 



Sitter (or de Sitter) space AdS2 (or dS2)- The phase space portrait of system (45) for 
the AdS2 case can be seen in Figure [13) 

5. Conclusions 

The gauge sector of the phase space geometry of a three dimensional volume preserving 
dynamical system is investigated. As a toy model we used the non dissipative part 
of the Lorenz system. Searching for possible valuable informations one may retrieve 
from this geometric picture, we found that a rich dynamical structure is revealed. 

This rich dynamical structure we are referring to, corresponds to different systems 
which in our application are the Duffing oscillator (Heisenberg case), the simple 
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Figure 13. The phase-space of system ( |13| > reduced on the hyperboloid Co for 
the case of AdS2- 



pendulum {SE{2) case), a uniformly charged rigid body in a uniform magnetic field 
{S0{3) case) and a mathematical 50(2, 1) formulation. All these turn out to be 
different representations of the same dynamics! Similar analysis can be performed 
to any three dimensional volume preserving dynamical system providing insight on a 
unified description of dynamics for specific systems. 

When the full dissipative system is considered the invariant manifolds of the 
conservative part can be used for the localization of the strange attractor. In 
addition, in section [TT] an elementary analogue of a 'spontaneous symmetry breaking 
mechanism' is proposed as an intuitive explanation for the formation of the attractor. 



For the S0{3) case, in section 4.3, the Lorenz system is physically identified with the 
rigid body system in a homogeneous magnetic field with external friction in every 
direction and a constant driven term in Z direction. A generalization of this rigid 
body system for arbitrary magnetic field's direction leads to a dynamical system with 
new strange attractors (figures [IT] and [l2| . 

Last, we note that the 5*0 (3) formulation could be useful for an angular 
momentum or spin-based quantization of the Lorenz attractor in accordance with 
the matrix formulation of ^ i9j. It may also be interesting to investigate the 
correspondence with the 5*0(2, 1) picture in this quantum realm. 
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Appendix A. 

Applying the Lagrange's multipliers method, we find the extremum of S{Xi^ Xs) 
under the constraint 5 = for: 



5 = (A - l)Xl + X| + (X3 - Ap)2 
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and (Xi,X2,X3) subject to the Lorenz flow Let call A = S. 

We are looking for a point {Xi^ Xs) and a Lagrange's multiplier a such, that: 
K = S{Xi^ X2^ Xs) is an extremum of S under the constraint A = 0. The constraint 
gives one equation: 

- a{X - l)Xf - Xl - hXl + 6(Ap)X3 = (A.l) 

Let extremize: f = S -\- aA. From: V/ = ^ V^* + a\/A = we get three more 
equations: 

Xi(A-l)(l-2acr) =0 (A.2) 
X2(1-2A) = (A.3) 



The system of Equations (A.1)-(A.4) has three distinct solutions: 

(i) {Xi,X2,X3) = (0,0,0) and a = 1/5 
which gives K = S{Xi,X2,X3) = {\pf 



(ii) (Xi,X2,X3) = (O,y'(Ap)20^,Ap2^) and a = 1/2 

which gives K — 5'(Xi, -^3) = (Ap)^ 4(6-1) • From X2 we see the constraint 
6 > 2. 



(iii) (Xi,X2,X3) = (^(Ap)2^^i|^,0, Ap^) and a = l/2a 

which gives K = X2, -^3) = — (Ap)^ 4cr(t-b) • ^^^^^ ^1 we see the constraint 

A<L 



Appendix B. 

Let a surface be defined by 5* = for some constant K and S = S{X\ If A = VS* • 
has the same sign at each and every point of S then is a localization surface. Let 
calculate A for Lorenz flow v ^ and 5' = (A - \)Xl + X| + (X3 - \pf". 

A = \JS-v^ (2(A - l)Xi, 2X2, 2(X3 - \p)) 

*(pX2 - aXi, -X1X3 + pXi - X2, X1X2 - 6X3) 
= - 2cr(A - \)Xl - 2X1 - + 26(Ap)X3 (B.l) 
We will prove that A has a constant sign for all repelling surfaces of Equation ( [31] ). 
• two sheeted hyperboloids A < 1 : S = —R^ = —{^p)'^ 4^^—^ 

Substituting: (A - 1)X^ = -X| - X| + 2(Ap)X3 - (Ap)^^^ in ^ we get: 

A = 2{a- 1)X| + 2(a - 6)X| - 2(Ap)X3(2a - 6) + 2a{Xpf ~ 

4<j((j — oj 

= 2('^-^)^^^ + 2(^ 

X {4X|((T - bf - 4(Ap)(2a - 6)X3((t -b) + {\pf{2a - bf) 

= 2{a - l)Xl + ^^-^ {2X3(a - 5) - (Ap)(2cT - 6)}^ > 
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V(Xi,X2,X3) G S and for a > 1, a > b. 
• cone A = : S = = {\pf = 



Substituting: Xf = X| + X| in ^BA\ we get: 

A = 2{a- 1)X| + 2(cr - 6)X| > 
V(Xi,X2,X3) G 5' and for a > 1, a > b. 



• parabolic surface A -oo : S = = (Ap)^ _j5^2 ^ 2^X3 



Substituting pX3 in (B.l) we get: 



A = (-Xi, 0, p) • (pX2 - aXi, -X1X3 + pXi - X2, X1X2 - 6X3) 

= aXf - bpXs = cjXI - ^-Xl = {a- ^-)X', > 
V(Xi,X2,X3) G 5' and for a > b/2. 

one sheeted hyperboloids A < : S = R^ = (Ap)^ 



Substituting: (A - l)Xf = -X| - X| + 2(Ap)X3 in (|b1]) we get: 

2(0- - 1)X| + 2(0- - 6)X| - 2(Ap)X3(2cr - 6) > 

V(Xi,X2) G 5 and X3 > and for cr > 1, a > b. Only the part X3 > of 
the surface is a repelling surface. From the parabolic surface we know that the 
attractor lies in X3 > 0. So these upper parts of one sheeted hyperboloids provide 
a lower boundary of the attractor. 
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